Two simple solutions 
of nonlinear Fokker - Planck equation 
for incompressible fluid 



Igor A. Tanski 
Moscow, Russia 
tanski. igor. arxiv@gmail. com 



Abstract 

In this article we derive two simple solutions of nonlinear Fokker - Planck equation for incompress- 
ible fluid and investigate their properties 
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1 Introduction 

In our previous work [1] we found following nonlinear Fokker - Planck equation for incompressible fluid: 

p = In dv\dv2dv^ = const, (1) 
Jv 

dn dn d(vjn) 1 dn dp d 2 n 
dt k dxk dvj pdvjdxj dvjdvj 

where 

n = n(t,xi,X2,x 3 ,vi,V2,v 3 ) - density; 

p - pressure; 

t - time variable; 

^1,^2,^3 - space coordinates; 
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vi,V2, V3 - velocities; 

a - coefficient of damping; 

k - coefficient of diffusion. 

In this work we constructed only stationary solution of the system (1-2) with zero average Maxwell 
velocities distribution and Pascal pressure field. In the present work we construct some another simple 
solutions, which generalize previous results. 

2 Zero pressure solutions of nonlinear Fokker - Planck equation for 
incompressible fluid 

First of all we consider zero pressure field 

p = 0. (3) 

In this case 

dn dn d(vjn) ^ d 2 n 
dt Vk dxk a dvj dvjdvj' 

This is simply standard Fokker - Planck equation, but solutions must satisfy additional equation (1). 
We know, that solution of equation (4) has following form (see [2]): 



+00 +00 +00 00 00 00 

n(t,xj,vj)= 12 12 12 12 12 exi p 



x (5) 



mi=— 00 m 2 =— 00 m 3 =— oo ni=0 n 2 =0 n 3 =0 

Xj4mim2m3?iiii2n3 < / , mim2m3niii2n3 > 

where 



tt ( m i v i\\ ( a i\ ( ( 4irimjk\\ 

^ mimmn2 n 3 = n exp (2m-( Xj - -) j exp (-- Vj ) H nj ^- ^ + —g— j j . (6) 

Comparing first term of (5) with constant p in the LHS of (1), we find: 

A mim2 m 3 000 = P l^^J • C* 7 ) 

(1) implies obviously, that for all another terms with m = n>2 = n 3 = coefficient is zero. 

A mim2 m 3 ooo = mi > or m 2 > or m 3 > 0; (8) 
For all another rtj we group terms with the same set of indices m; and equal sums of indices nf. 
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G s (mi,m 2 ,m 3 ) 



J2 exp 

ni+ri2+n3=s 



<(•£**' ©',§(?)') 



xA 



mim2m3nin2n3V^mirri2m3nin2«3 ) 



where group index 



(9) 



E 



n-j- 



(10) 



We considered such groups in our works [3] (s = 0) and [4] (s = 1). 

For each group G s amplitude decreases with time as exp(— ast). Therefore we can carry out time 
dependent multiplier before the sum: 



G s (mi, 7712,777,3) = exp 



x ^ ' j4 mim2m3 ij in2fl;s (/) mim2m3 ij in2n3 , 

n-i+n,2+n-3=s 



Let us denote the integral: 



m = exp (-2iri 



exp -— v 2 H, 



a-i a l \ 2k 3 1 



Vi + 



4nimjk 



2/7 - 



Equation (1) gives following orthogonality condition for each group: 

n 1 +n 2 +n 3 =s 

To calculate the integral (1) we perform following substitutions: 



(11) 



(12) 



(13) 



a / 2 <iirimkv\ a 
2k \ + a 2 a J = ~2k 



v + 



2mmk\ 2 /2nmk\' 

-) + v^r) 



a* a 



I a ( 2nimk 

i = HTk\ v + ^T 



This gives: 



(14) 
(15) 



■ exp 



a / 2"7T777,/c \ 2 

2k V a 2 a ) 



exp 



, , , a 2-imk \ 

< 2 H n [£ + J — —2—)dZ. 

2k a z a J 



(16) 



The addition theorem for Hermit polynomials is (see [5]): 
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n 

' n 



H n ( Zl + z 2 ) = 2"™/ 2 1 H k (ziV2)H n - k (z 2 V2). (17) 



Orthogonality condition for Hermit polynomials is: 



r exp (-x 2 ) H m (x)H n (x)dx = 5 mn ^2 n n\. (18) 

Apply (17) to (16) 

«• (^") ^ 2 -" /2 £ (I) *«^«~ (^) ^ < 19 > 



= ^.(^ (^) + 2 -„ M *„-, (^) + .... 

In the first version of this article we applied (18) to (19) an got erroneous result. In fact we can not 
apply (18) to (19), because argument of H n in (19) is (£V2), not (£) - as in (18). So addition theorem in 
the form (17) is useless for integral I mn calculation. 

Instead of this there are three correct ways to obtain integral I m n value: 

- direct use of generating function of Hermite polynomials. We consider this way in APPENDIX 1. 

- another form of addition theorem. We consider this way in APPENDIX 2. 

- use of Fourier transform. We consider this way in APPENDIX 3. 

This last way we used in our work [2]. This makes our error in the first version of this paper even 
more inexcusable. 



The result is: 



2irk (2k\ n l 2 
Imn = \ exp 



a 



k /2vrm\ 2 
2a \ aa ) 



/2 ™)\ (20) 



V aa 



Carry out all multipliers, which depend only on m; or s = J2 n i- This gives following orthogonality 
condition for the group: 

\- A ( m A ni f m A n2 f m A n3 -n fou 

7 j /i -m 1 m 2 m 3 n 1 n 2 n 3 — U S > U. (Zl) 

ni+n2+n , i=s \aiJ \a 2 J \a 3 J 

For example for the group s = 1 we have following condition: 

777/1 TFl2 

^4mim2m3l00 l~ ^-rnirr^rn-iOW I" A mim2m3 o@i =0 S = 1. (22) 

cti a 2 

This condition means physically, that divergence of the average velocity field is zero. When initial 
velocity field satisfy the set of such conditions, it can be prolonged for all times tO with zero pressure 
field. 
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3 Average velocities 



It is interesting to derive expressions for average velocities. For this purpose we denote: 




In the first version of this paper we tried to calculate J mn value using addition theorem (17). But, 
as before, (17) is inconvenient for this purpose. The correct way of J mn value deduction is given in 
APPENDIX 4. 

So we have 



Jmjnj 



2irk f2k\ n ^ 2 



a \ a J 



exp 



k I 2nnij 
2a \ aa,j 



1 /r, • \ n 3+ 1 /r, • \ n J~ 1 

- ' - ' 1 + nj 1 '- 1 



a \ 



aciA 



(24) 



2irk 

JmjO = ~y — eX P 



k I 2mn~ 



2a \ era,- 



k / 2mm j 
a \ acij 



(25) 



In this way we get final expressions for pui 



pu\ = exp 



3 /0 ^2 3 

-t a n j + k 



■j 



a I " \ cij 

j=i \ j 



3=3 

TT I ^ - m i 

j _[ exp 2m — -Xj 

i=i 



(26) 



pu2 = exp 



J=3 



TT / m 7 

1 1 exp 2m — - 



X,' x 



a. 



x £ A 

n 1 +n 2 +n 3 =s 



(27) 



pu 3 = exp 



t la "£n j + k (^) E 



\ 3=1 

x 2 A 

ni+n2+"3=s 



V a J . \ 

.7 = 1 



CLj 



Q exp ( 27ri — =^-Xj 



n 1 Yraimznxnino, Im\n\ IrnzniJmznz ■ 



(28) 



In the next section we consider solution 



Anoioo — 1; Anooio — 0; Anoooi ~~ ®- 



(29) 



In this case 
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n = p 



{ r ( 27r \ : 

+ exp —t [a + k 

\ \aa 2 J 

Coefficients are equal to 



a \ 3 / 2 
2^k) 6XP 



a 

2k 1 



exp l (x 2 

\ a 2 



a 



exp 



a 



2^ 



(30) 



^2fc r 2-Kk 

Joi = Vtt — , ho = \ exp 

a \ a 



a ( 2-irk \ 2 



2k \a 2 a 2 



, ho — 



l 2ixk 



a 



1 {2irhy 
pui = —j=\ exp 



a ( 2-nk \ 2 



2k \a 2 a 2 



exp 



-t \a + k (—} 

\ \aa 2 J 



f2m 

exp x 2 

\ a 2 



U2 = u 3 = 0. 



(31) 

(32) 
(33) 



This is time dependent shift along X-axis with amplitude depending on Y space coordinate. Note, 
that u\ is periodic on x 2 and average u\ on x 2 is zero. 

4 Constant pressure gradient solutions of nonlinear Fokker - Planck 
equation for incompressible fluid 

The case of constant pressure gradient is also very simple. Let us denote 



dp dp dp 

— = ah; — = 0; — = 0. 

OXi ox 2 OX3 



Equation (2) take the form 



dn dn divkn) dn 
+ Vk~ a — h— — = k 



d 2 



■71 



dt dxk dvk dvi dvkdvk 

Perform in (35) the substitution 

v[ = v± — h; vi = v[ + h; x\ = x\ — ht; x\ = x\ + ht. 

and (35) transforms to zero pressure case (4). 

This follows also from eq. (68) of [1]. More specifically, we use the symmetry group 

6 = («/i-/i>| + /!w|; + AW|i 

where 



(34) 



(35) 



(36) 



(37) 



/i=t; /i = l; /C = 0, 



(38) 
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that is 



d d d 
vz = ax— + — +t—. 
op ou ox 



(39) 



The invariant solutions are 



n(t, x, y, z, u, v, w) = no(t, x + ht, y,z,u + h, v, w). 
where no is solution for zero pressure case. 



(40) 



p(t, x, y, z) = ah x. 



(41) 



For the solution (29-30), which does not depend on X coordinate, we get in this way Poiseuille type 
solution. 

In this case 



n = exp 



. 2m v 2 , 
exp (x 2 ) I x 

a>2 a 



-t (a + k f— \ 
\ \aa 2 J 



x exp ( -— ((«! + h) z + vi + v£) \ ^ ( ^ — ( Vl + h) ) ; 



(42) 



1 /2irk\ 2 

pu\ = —ph H — exp 

\ a J 



7T 



a / 2nk 



2k \a 2 a2 



exp 



-t (a + k (—\ 
\ \aa 2 J 



2iri 
exp l x 2 

a.2 



X2) ■ (43) 



U2 = U3 = 0. 



(44) 



We see, that constant increasing of p along X axis results in adding constant velocity in the opposite 
direction - quite as for Poiseuille flow. 



DISCUSSION 



In this article we derive two simple solutions of nonlinear Fokker - Planck equation for incompressible 
fluid and investigate their properties. They are: flows with zero pressure and flows with constant pressure 
gradient. This short list of solution will be useful for further investigations. 
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APPENDIX 1 



We start from generating function of Hermite polynomials (see for example [5]): 

00 t n / \ 

H n (z) — ! = exp [-t 2 + 2ztj . 



n=0 



(API - 1) 



Substitute z = z\ + z 2 



t n 

V H n (zi + z 2 ) — \ = exp ( -t 2 + 2(zi + z 2 )t 



n=0 



Multiply both sides by e Zl and integrate: 



e" z i H n( z i + z 2 ) — \dzi = / e _2 i exp (-t 2 + 2( Zl + z 2 )t) dz\ 

-co n=0 n J -co ^ ' 

exp(2z 2 t) f°° e-^-^dz! = ^Fexp(2z 2 t) = 



n=0 



(we used Poisson's integral value). Compare coefficients by t n and get 



(API - 2) 



(API - 3) 



/ e z iH n (zi + z 2 )dzi = V^(2z 2 ) 

J —CO 



(API - 4) 



To get value of I mn integral we substitute z 2 = a/^ 2 "^ - We have by definition 



L mn — \/ exp 

a 



a /'2irmk\ 2 



2k V a 2 a 



exp 



a 2 r Kimk 
2k a 2 a 



(API - 5) 
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According (AP1-4) 



'27T/C 



a 



■ cxp 



a f 2irmk 
2k V a 2 a 



2a\ n / 2 (2iximk\ n 
T 



a 2 a 



So we have 



k , (2k\ n l 2 



— exp 

a 



k /27rm\ : 
2a V aa / 



k /2vrm\ 2 
2a V aa / 



27rmV 
aa / 



/ 2mm \ 1 
\ aa J 



(API - 6) 



(API - 7) 
(API - 8) 



APPENDIX 2 



Use once again (AP1-2). This time we multiply both sides by e z ^H m (zi) and integrate: 



e - 2 i# m (zi) ^ fl- n (zi + z 2 ) — \dzi = / e~ 2 i# m (zi)exp (-t 2 + 2(zi + z 2 )t 

-OO n=0 ^ °° 

= exp (2z 2 t) / i7 m (zi)e _(2:i_t)2 dzi = exp (2z 2 t) / F m (x + t)e~ x2 dx 

J —CO J —CO 



dzi = (AP2 - 1) 



— OO 



exp (2z 2 t) v^F (2t) m = v^F (2i) m f] ^ 



fc=0 



Compare coemcients by 



/*°° _ t 2 , . , J n < m 
J-oo 6 Hm[t)Hn[t + X)M - { ^2 n x n - m n\/(n -m)\ n > m ' 


(AP2 


-2) 


This integral gives coemcients of following Fourier decomposition: 


n / \ 

#„(* + *) = £(?) (2x) n - k H k (t) 
k=0 v ' 


(AP2 


-3) 



This is another form of addition theorem (compare with (17)), which we can use to calculate /„ 
integral value. 

Another way to deduct (AP2-3), which considers (AP2-3) as Taylor series, is given in [6]. 
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APPENDIX 3 



By definition 

According to [2] Fourier transform of (f) mn is 



ran — / Ymn 

oo 



dv. 



(AP3 - 1) 



k ( 2-Km\ 
za V aa / 



27rm\' 



aa 



= F mn (q). (APS - 2) 



According to definition of Fourier transformation we could express the value of integral of any function 
through the value of its transform in zero point: 



(aj GXP 



27rF m „(0) = 2ttx 

k ( 2vrm\ 2 ' 
2a V aa / 



2vrmY 



aa 



a = 0. 



This gives 



(AP3 - 3) 



or 



fc /2vrm\ 2 
2a V aa / 



2nm\ 1 



aa 



- 4) 



2tt£; /2fc\ n / 2 



a \ a 



exp 



fc /2vrm\ 2 
2a \ aa / 



(2mm\ n 
\ aa J 



(AP3 - 5) 



APPENDIX 4 



By definition 



,v dv. 



(APA - 1) 



According to definition of Fourier transformation we could express the value of first moment of any 
function through the first derivative value of its transform in zero point: 



J mn = 2iriF* mn (0) = 2TTiF mn (0) 



k ( 2irm\ n 

— U+ + 7 n \~ 

a V aa J (g_ j 



(APA - 2) 



q=0 
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We have finally: 



Jn 



2irk /2k\ n/2 



— exp 

a 



k /2vrm\ 2 
2a \ aa J 



k /2irim\ n+1 



a \ aa 



+ n 



( 2itim \ 
V aa ) 



n-l 



(APA - 3) 
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